Entanglement in a hardcore-boson Hubbard model 
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The entanglement in a Hubbard chain of hardcore bosons is investigated. The analytic expression 
of the global entanglement in ground state is derived. The divergence of the derivative of the global 
entanglement shows the quantum criticality of the ground state. For the thermal equilibrium state, 
the bipartite and the multipartite entanglement are evaluated. The entanglement decreases to zero 
at a certain temperature. The thermal entanglement is rapidly decreasing with the increase of the 
number of sites in the lattice. The bipartite thermal entanglement approaches a constant value at 
a certain number of sites while the multipartite entanglement eventually vanishes. 
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I. INTRODUCTION 



Quantum entanglement is an important property 
which plays an essential role in the quantum informa- 
tion processing [3, HI, 0]. There are different types of 
entanglement, such as bipartite, multipartite, block 
and localizable entanglement, etc Q. The relative en- 
tropy of entanglement [|| and the entanglement of for- 
mation Q can be used for qubits. Arbitrary bipar- 
tite entanglement can be assessed by the "negativity" 
L 8|. In recent years, the entanglement in Heisenberg 
models of finite s yste ms of spins has been investigated 
EJMEI EE EEIIE EE EE G3]. The anisotropy effect 
ldlllj , multi-dimensional lattices [HJ , lon g -ran ge in ter- 
actions [l3j], multipartite entanglement 0, EEUE E3 
have all been studied in Heisenberg models. The en- 
tanglement in solids can also be witnessed by the mag- 
netization [TU and the thermal energy [IE HE HI| . In 
solids, there is a characteristic temperature below which 
a thermal entangled state can be obtained. The ef- 
fects of quantum entanglement have been detected in 
the experimental situation (22J. Some attention has 
been drawn to the entanglement in the models of infi- 
nite spin systems [23l . |24| . Owing to the quantum nonlo- 
cal correlations, the connection of the entanglement and 
quantum criticality i n sp in systems has been discussed 
[H HE US HE HE M HE HE S [p . The quantum crit- 
icality can be shown by the entanglement of the ground 
state JEHU- Next to Heisenberg models, boson Hubbard 
models have been extensively used to study the metal- 
insulator transition [35] . The Hubbard chain of hardcore 
bosons is one of the simplest models which can embody 
such quantum criticality. It is of interest to investigate 
the entanglement in the hardcore boson Hubbard chain. 

In this paper, the entanglement of a hardcore boson 
Hubbard chain is studied for the ground and thermal 
states. In Sec. II, the global entanglement measure of 
the ground state is analytically expressed by the magne- 
tization. The critical properties of the transition between 



the Mott-insulating phase and the superfluid phase are 
shown by the derivative of the global entanglement at 
ground state. In Sec. Ill, the "negativity" of any two 
sites in the chain is derived. The multipartite entan- 
glement at finite temperature is calculated. The effects 
of the number of lattice sites, the temperature, and the 
chemical potential are investigated. A discussion con- 
cludes the paper. 



II. GLOBAL ENTANGLEMENT OF GROUND 
STATE 

Recently, the quantum entanglement of strongly corre- 
lated spin systems has been extensively studied. One of 
the interesting focuses was the relation of quantum phase 
transitions and entanglement in the ground state. As is 
well known, the phenomenon of quantum phase transi- 
tions describes the global properties of the ground state. 
Therefore, the measurement of global quantum entan- 
glement is relevant to the investigation of the quantum 
phase transition. For a typical case, the spin version of 
the one-dimensional boson Hubbard is expressed by 
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where a\ and asj are the creation and annihilation oper- 
ator for bosons on the site i of the lattice, n a i = a\a.i 
is the number operator. The parameter w > allows 
hopping of bosons from one site to another, /i represents 
the chemical potential of the bosons, and U denotes the 
possible repulsive interaction among bosons on each site. 
For simplicity, only the nearest-neighbor pairs (ij) are 
considered. The off-site and long-range repulsive inter- 
actions are neglected. If the repulsion is very strong, 
when U — > oo, there is only one boson at each site. The 
model is reduced to the simplest one of a hardcore boson 
Hubbard chain, which can also be written as a magnetic 
model of S = | spin with pairwise interaction. The rela- 
tion of = 1— 2a\a,i is satisfied. In this following, the hard- 
core boson Hubbard chain is studied. It is known that 
this model is equivalent to the spin- 1/2 Heisenberg XX 
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chain with the ferromagnetic interaction —w/2 and the 
external magnetic field p. Through the Jordan- Wigner 
and Fourier transformations, the Hamiltonian H can be 
exactly expressed by 



global entanglement E — 0. It is found that the quantum 
criticality in the hardcore boson Hubbard chain can be 
shown by the global entanglement at the ground state. 



H 



zZ ekc l ck 

k 



(2) 



III. ENTANGLEMENT AT FINITE 
TEMPERATURE 



-2wcos2j^ 



p, 



L is the number of sites 
are the Fourier-transformed 



where 

on the lattice, and c ] k ,Ck 
fermionic operators. 

To show the relation of entanglement and quantum 
criticality, the global entanglement E of the ground state 
can be introduced by [3(| [33] , 



E = 2[1 
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(3) 



where tr(p 2 ) is the trace of the reduced density matrix p t 
on the ith site of the ground state. Here, the eigenstates 
of of are assumed to be {|0)<, |1)<}. The reduced density 
matrix is given by pi — where M=j- tr(pi<jf). 

Therefore, the global entanglement E of the ground state 



can be just obtained by 



E = 1 — M 



For the values of \p\ < 2w, there is partial occupation 

of sites at the ground state and M — — cos -1 ^ 1. 

When \p\ > 2w, \M\ = 1. The global entanglement 
of the ground state is plotted as functions of the hop- 
ping coefficient w and the potential p in Fig. 1(a) when 
the number of sites is L = 10 4 . It is shown that the 
global entanglement E exists at the ground state if the 
potential satisfies \p\ < 2w. The values of E are decreas- 
ing in \p\ and then drop to zero when \p\ — 2w. When 
p —> 0, it is the maximally entangled ground state which 
is exactly the Greenberger-Horne-Zeilinger state of the 
form \i>) g = -^=(|010 - - - 1> + |101 •■•())). If the values 
p < —2w or p > 2w, the ground state is an unentan- 
gled pure state 00 • • ■ 0) or 11 • • ■ 1). According to [35j ] . 
there are two different kinds of phases namely the Mott- 
insulating phase and the superfluid one. It is found that 
the transitions between them occur under the condition 
p = ±2w. The global entanglement always exists in the 
superfluid phase while there is no entanglement in the 
Mott-insulating phase. To clearly demonstrate the phe- 
nomenon of a quantum phase transition in the ground 



state, the derivative of the global entanglement E u 
tained 



is ob- 



dE 
dp 



TT 2 ^4w 2 



pcos- 1 -^ 

1 2w 



I/'- 



< 2w. 



the divergence of ^ . The quantum criticality is depicted 
in Fig. 1(b) when the hopping coefficient is chosen to 
be w = 1. If p — > 2w + , the derivative |p = for the 



When the potential p — ► 2W-, ^ 



(5) 

— co which reveals 



The thermal equilibrium state is p(T) = e~ H/>kT /Z, 
where Z is the partition function at finite temperature 
T and k is the Boltzmann constant. For convenience, 
both the Boltzmann constant k and the Planck constant 
h are assumed to be one. Because the Hubbard chain of 
hardcore bosons is equivalent to the spin- 1/2 Heisenberg 
XX chain, the reduced density matrix pij on any two 
sites i and j can be expressed by the correlation function 
K?j* = tr{pcrfcr^) : (a = x, y, z). In the Hilbert space of 
{|00)i 3 , |01)y, |10)y, |ll)jj}, the expression for p^ can be 
obtained 



Pij 
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where u = \{K$ + 2M + l),v= \{K?f -2M+\),w = 
(4) |(1 - Kff), and t = \Kf?. The analytical calculations 
of correlation functions K°j a and the magnetization M 
are straightforward given. For the number of sites L, the 
magnetization is given by M = — X)g=i tanh(e g /2T). 
The two-site correlations can be by 



T\ xx 

ij 



Gi Gq 
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(7) 



where the parameter r=\j — i\ 
tion distance between two sites, 



is the separa- 
and the item 



G r =G_ r = i Y^=i cos(2-Kqr/L) tanh(e g /2T). 

The thermal entanglement in the chain can be inves- 
tigated by the negativity N. Based on the separability 
principle, the negativity N can be used to quantify the 
bipartite entanglement between two sites The nega- 
tivity N is introduced by 



N(p) 



(8) 



where Xi is the ith negative eigenvalue of p T which is the 
partial transpose of the mixed state p. From the sepa- 
rability of quantum states, the partial transpose matrix 
p T has nonnegative eigenvalues if the states are unentan- 
gled and the value N(p) — 0. In the Hubbard model of 
hardcore bosons, the negativity of the two-site entangled 
state is given by 



N(pij) = \u + v- {u - v) 2 + At 2 \ 



(9) 
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Thus, the entanglement on any two sites can be calcu- 
lated numerically through Eqs. ([T])-©- The negativity 
N is plotted as a function of site number L in Fig. 2(a) 
when the temperature is T = 0.5, the hopping coefficient 
w = 1, and the potential fi = 0.2. It is found that the 
pairwisc entanglement always exists whatever the num- 
ber L of sites is. There is no thermal entanglement on 
two sites when the separation distance is r > 2. The 
value of N is rapidly decreasing in L, and then reaches a 
constant value at L = 16. This result illustrates that the 
thermal entanglement can be detected in real solids of a 
very large number of particles. 

The multipartite entanglement El for the thermal 
states of systems with an even number of sites L can 
be introduced 

2 L -1 

E L (p) = max{0, v a - ^ Vj } (10) 

3=1 

where {Vj}?—^ 1 is the spectrum of the operator 

y/ pUpU -1 in decreasing order, U is an anti-unitary 
time reversal operator and can be written as U = 
[TljL 1 (— i<Tj)]T, t is the complex conjugate operator. 
The multipartite entanglement El is plotted in Fig. 2(b). 
In Fig. 2(b), the value of El is decreased with the in- 
crease of the number of sites L. The multipartite entan- 
glement El vanishes at L = 10. However, the bipartite 
entanglement still exists in this case. 

It is also very interesting to study the effects of the 
temperature T and the chemical potential /i on the ther- 
mal entanglement. By the analytical expression for the 
bipartite entanglement in Eq. ([9]), the negativity N is 
plotted as a function of T and /i in Fig. 3(a) for w = 1 
and L = 10 4 . When the temperature T and the chemical 
potential \fi\ are increased, the negativity N decreases. It 
is found that the values of the negativity N are symmetric 
about the chemical potential p. When /j, — ► 0, the value 
of N is maximal. For a definite chemical potential, the 
bipartite entanglement disappears at a certain tempera- 
ture T c . It is clear that the values of T c can be increased 
by decreasing the chemical potential \p\. It is seen that 
the entanglement can be detected at low temperatures in 



solids. The multipartite entanglement El for the ther- 
mal state of L = 6 is plotted in Fig. 3(b) when w = 1. It 
is shown that the values of El are also symmetric about 
the chemical potential p and decreased with the increase 
of \p\. The ground state for \/j,\ < 0.5 is just the max- 
imally entangled GHZ state -^(|010 • • • 1} + 1 101 ■ • -0}). 
The values of El are declined with the temperatures T 
and vanishes at about T = 0.6. 



IV. DISCUSSION 

The entanglement in a hardcore boson Hubbard chain 
at ground and thermal equilibrium states is investigated. 
The global entanglement at ground state is analytically 
expressed by the magnetization. When the potential 
fi — ► 0, the maximally entangled Greenberger-Hornc- 
Zeilinger state can be obtained. The quantum critical- 
ity is revealed by the divergence of the derivative of the 
global entanglement at ground state. In the parame- 
ter plane of Mott-insulator and superfluidity phases, it 
is found that the entanglement exists in the SF phase 
while there is no entanglement in the MI phase. The bi- 
partite entanglement between any two sites is deduced 
by the negativity. For a very large number of sites L, 
the pairwise entanglement can always exist. When the 
number L increases, the negativity decreases rapidly and 
then reaches a constant value at a certain number of site. 
While the multipartite entanglement will decrease to zero 
with increasing L. The thermal entanglement vanishes at 
a certain temperature and is decreased with the increase 
of the potential p. It is shown that the entanglement can 
be detected at low temperature in real solids of a large 
number of sites. 
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Figure Captions 
Fig. 1 

(a) The global entanglement E of the ground state is 
plotted as a function of the hopping coefficient w and 
the potential \i. The number of sites is L — 10 4 ; (b)Thc 
derivative ^ is plotted to show the quantum criticality 
of the ground state. 

Fig. 2 

The thermal entanglement is plotted as a function of 
the number of sites when w = l,fj, = 0.2, and T = 0.5. 
(a) The pairwise entanglement of the negativity N; (b) 
The multipartite entanglement El. 

Fig. 3 

The thermal entanglement is plotted as a function of 
the potential /x and the temperature T when the hopping 
coefficient is w = 1. (a) The pairwise entanglement for 
L = 10 4 ; (b) The multipartite entanglement El for L = 
6. 
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